Dynamics of pion production in heavy-ion collisions around lA GeV energies 
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Within the framework of the improved isospin dependent quantum molecular dynamics (ImlQMD) 
model, the dynamics of pion emission in heavy-ion collisions in the region of 1 A GeV energies as a 
probe of nuclear symmetry energy at supra-saturation densities is investigated systematically. The 
total pion multiplicities and the -k~ /tt'^ yields are calculated for selected Skyrme parameters SkP, 
SLy6, Ska and SIII, and also for the cases of different stiffness of symmetry energy with the parameter 
SLy6. The influence of Coulomb potential, symmetry energy and in-medium pion potential on the 
pion production is investigated and compared each other by analyzing the distributions of transverse 
momentum and longitudinal rapidity and also the excitation functions of the total pion and the 
tt^/tt^ ratio. The directed flow, elliptic flow and polar angle distributions are calculated at the 
cases of different collision centralities and also the various stiffness of the symmetry energies. A 
comparison of the calculations with the available experimental data is performed. 
PACS number(s): 25.75.-q, 13.75.Gx, 25.80. Ls 



I. INTRODUCTION 

Heavy-ion collisions induced by radioactive beam at 
intermediate energies play a significant role to extract 
the information of nuclear equation of state (EoS) of 
isospin asymmetric nuclear matter under extreme con- 
ditions, such as high density, high temperature and large 
isospin asymmetry etc. Besides nucleonic observables 
such as rapidity distribution and flow of free nucleons or 
light clusters (such as deuteron, triton and alpha etc.), 
also mesons emitted from the reaction zone can be probes 
of the hot and dense nuclear matter. The energy per nu- 
cleon in the isospin asymmetric nuclear matter is usually 
expressed as E{p, S) = E{p, S ^ 0) + Esyin{p)S^ + 0{6^) 
in terms of baryon density p ^ pn + Pp, relative neutron 
excess 5 = {pn — Pp)/(,Pn + Pp), energy per nucleon in a 
symmetric nuclear matter E[p,5 = 0) and bulk nuclear 

symmetry energy ^;sym = \ l<5=o- In general, two 

different forms have been predicted by some microscopi- 
cal or phenomenological many-body approaches. One is 
the symmetry energy increases monotonically with the 
baryon density, and the other is the symmetry energy 
increases initially up to a supra-saturation density and 
then decreases at higher densities. The difference of the 
symmetry energy at supra-saturation densities predicted 
by transport models is huge. It is not only in understand- 
ing the reaction dynamics, the high-density behavior of 
the symmetry energy also has an important application 
in astrophysics, such as the structure of neutron star, the 
cooling of protoneutron stars, the nucleosynthesis during 
supernova explosion of massive stars etc [l|. With the 
establishment of high-energy radioactive beam facilities 
in the world, such as the CSR (IMP in Lanzhou, China) , 
FAIR (GSI in Darmstadt, Germany), RIKEN (Japan), 
SPIRAL2 (GANIL in Caen, France) and FRIB (MSU, 
USA) , the high-density behavior of the symmetry en- 
ergy can be studied more detail experimentally in the 
near future. 



Based on recent analysis of experimental data associ- 
ated with transport models, a symmetry energy of the 
form Esymip) ~ 31.6(p/po)^ MeV with 7 = 0.69 - 1.05 
was extracted for densities between O.lpo and 1.2pQ 
Theoretically, the symmetry energy at supra-saturation 
densities can be investigated by analyzing isospin sen- 
sitive observables that are emitted in the high-density 
region, such as the neutron/proton ratio of emitted nu- 
cleons at mid-rapidity, 7r~/7r+, S]^/E+ and K°/K~^ ra- 
tios etc Q. Heavy-ion collisions around lA GeV ener- 
gies can reach a 2-3po nuclear matter. The emission of 
pion in heavy-ion collisions in the region 1 A GeV is 
especially sensitive as a probe of symmetry energy at 
supra-saturation densities. Recently, a very soft symme- 
try energy at supra-saturation densities was pointed out 
through fitting the FOPI data [i| by using IBUU04 model 
The inverse results were also reported by the Catania 
group with the RBUU model, in which a hard symmetry 
energy in the high-density region resulted in the large 
values of the 7r~/7r"'" and K^/K^ ratios 0,13]. Further 
investigations of the pion emissions in the 1 A GeV region 
are still necessary by improving transport models or de- 
veloping some new approaches. The ImlQMD model has 
been successfully applied to treat dynamics in heavy-ion 
fusion reactions near Coulomb barrier and also to de- 
scribe the capture of two heavy colliding nuclides [1, 0| . 
Some further improvements of the ImlQMD model have 
been performed in order to investigate the pion dynamics 
in heavy-ion collisions [13, U | . 



A systematic comparison of the ImlQMD results and 
the available experimental data is performed in this work. 
Calculations are focused on the pion production in heavy- 
ion collisions around lA GeV energies. Some observables 
sensitive to the high-density behavior of the symmetry 
energy are investigated and discussed. The paper is or- 
ganized as follows. In Sec. II we give a detailed de- 
scription of the ImlQMD model. Calculated results of 
pion dynamics in heavy-ion collisions and constraining 



2 



the symmetry energy at supra-saturation densities are 
given in Sec. III. In Sec. IV conclusions are discussed. 



consists of the relativistic energy, the effective interaction 
potential energy and the momentum dependent part as 
follows: 



II. MODEL DESCRIPTION 



The same as the QMD ^ or IQMD model the 
wave function for each nucleon in ImlQMD is represented 
by a Gaussian wave packet 



(27rL)3/4 



exp 



X exp 



ipi{t) ■ r 



(1) 



Here ri{t), Pi{t) are the centers of the ith nucleon in the 
coordinate and momentum space, respectively. The L 
is the square of the Gaussian wave packet width, which 
depends on the mass number of nucleus. The total N- 
body wave function is assumed as the direct product of 
the coherent states, where the anti-symmetrization is ne- 
glected. After performing Wigner transformation for Eq. 
(1), we get the Wigner density as 



(2) 



(7) 



Here the Pj and rrii represent the momentum and the 
mass of the baryons. 

The effective interaction potential is composed of the 
Coulomb interaction and the local interaction 

Uint = Ucoul + Uloc- (8) 

The Coulomb interaction potential is calculated by 

Ucoui - i E ^erfin,/VlL) (9) 

where the is the charged number including protons 
and charged resonances. The r^- = jr^ — | is the relative 
distance of two charged particles. 

The local interaction potential energy is derived di- 
rectly from the Skyrme energy-density functional and ex- 
pressed as 



with 



■ exp 



(r-r,(t))2 (p-p,(t))2.2L 
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The density distributions in coordinate and momentum 
space are given by 



=y /(r,P,i)dp 



(27rL)3/: 
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The local potential energy-density functional reads 
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(5) 



respectively, where the sum runs over all nucleons in the 
reaction systems. 

The time evolutions of the baryons and pious in the 
system under the self-consistently generated mean-field 
are governed by Hamilton's equations of motion, which 
read as 

dH . dH 

Pi = - — ' = — ■ (6) 



dpi 



Here we omit the shell correction part in the Hamiltonian 
H as described in Ref. f3| . The Hamiltonian of baryons 



where the pm Pp and p — Pn + Pp are the neutron, 
proton and total densities, respectively, and the 5 — 
{Pn ~ Pp)/{Pn + Pp) is the isospin asymmetry. Here, all 
the terms in the Skyrme energy functional are included 
in the model besides the spin-orbit coupling. The coef- 
ficients a, j3, 7, gsur, ffsur-i 9t are related tojfche Skyrme 
parameters 7 ^11*2,^3 and Xq,Xi,X2,x^ as IQI], 



1-3 7 



(12) 
(13) 



a _ 3 fi 
2 "8^°^°' TT"7"16''"' 

^TT" = ■^(9*1 ^ - 4X2^2)^0, 

1 d4 

^ = -^[3^1(2X1 + 1) + t2{2x2 + l)]po, (14) 

5r = ^ QtT^^ ' (3ti + hh + ^X2t2)pl'^. (15) 
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The parameters of the potential part in the bulk sym- 
metry energy term are also derived directly from Skyrme 
energy-density parameters as 



yrn 

1 

2/3 



48 



-sym 



^0^') pT[Shx,^t2i5x,+4)]. (16) 

The momentum dependent term in the Hamiltonian is 
taken as the same form in Ref. [15] and expressed as 



U„ 



1 Y: ^[ln(.(p,-p,.f + 

2 po 



with 



Pij 
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exp 
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(17) 



(18) 



which does not distinguish between protons and neu- 
trons. The parameters 5 and e were determined by fitting 
the real part of the proton-nucleus optical potential as a 
function of incident energy from the experimental elastic 
scattering data. The last two terms in Eq. (11) origi- 
nate from the momentum-dependent part of the Skyrme 
interaction which, in principle, should contribute to the 
nucleon effective mass [16, 17]. However, in the ImlQMD 
model their contribution is taken into account in the 
potential energy-density functional only and not in the 
momentum-dependent interaction given by the term of 
Eq. (17) solely. Therefore, the effective (Landau) mass 

in the model is fixed as ml^ = + with 

the free mass m at Fermi momentum p = p^ [18], which 
produces the value of 7n*^/m = 0.77 for all sets of the 
model parameters. 

In Table 1 we list the ImlQMD parameters related to 
several typical Skyrme forces after including the momen- 
tum dependent interaction. The parameters a, /3 and 
7 are redetermined in order to reproduce the binding 
energy (£'^=-16 MeV) of symmetric nuclear matter at 
saturation density po a-nd to satisfy the saturation prop- 
erties by relation ^^Jj^ lp=pQ=0 for a given incompress- 
ibility. Combined Eq. (16) with the kinetic energy part, 
the symmetry energy per nucleon in the ImlQMD model 
is given by 

2/3 

^sym (p) 



'32^ 



-TT p 
P_ 

Po 



Po 



5/3 



(19) 



More clearly compared with other transport models, such 
as IBUU04, RBUU, and IQMD etc, the bulk symmetry 
energy can be expressed as 



E. 



sym 



2/3 



pot 



(20) 




FIG. 1. The density dependence of the nuclear symmetry en- 
ergy for different Skyrme forces SkP, Sly6, Ska and SIII (su- 
persoft), and for the cases of soft, linear, hard and superhard 
trends. 



The value 7s = 1 is used in the IQMD model ^S^- In 
Fig. 1 we show a comparison of the nuclear symmetry 
energy at the situations of different Skyrme forces SkP, 
SLy6, Ska and SIII calculated by Eq. (19), and also the 
cases in Eq. (20) of ^s— 0.5 (soft), 1 (linear), 2 (hard), 
3 (superhard) with the coefficient Cpot =38 MeV, which 
have the value of the symmetry energy i?sym(po)=31.5 
MeV at the saturation nuclear density po=0.165 fm~'^. 
The parameter SLy6 is used in the calculations by de- 
fault. 

The pion is created by the decay of the resonances 
A(1232) and N*(144G) which are produced in inelastic 
NN scattering. The reaction channels are given as fol- 
lows: 

iViV^TVA, NN^NN\ NN AA, 

A o Nn, N* o Nn, NN iViV7r(s-state). (21) 

The cross sections of each channel to produce resonances 
are parameterized by fitting the data calculated with the 
one-boson exchange model fl^. In the 1 A GeV region, 
there are mostly A resonances which disintegrate into a 
TT and a nucleon, however, the N* yet gives considerable 
contribution to the high energetic pion yield. The energy 
and momentum dependent decay width is used in the 
calculation [lOj]. 

Analogously to baryons, the evolution of pions is also 
determined by the Hamiltonian, which is given by 



where the Pj and represent the momentum and the 



E ( + \/p' + < + ReK'^*(p„ p.) ]m) 
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TABLE I. ImlQMD parameters and properties of symmetric nuclear matter for Skyrme effective interactions after the inclusion 
of the momentum dependent interaction with parameters 5=1.57 MeV and e=500 c^/GeV^. 



Parameters 


SkM* 


Ska 


SIII 


SVI 


SkP 


RATP 


SLy6 


a (MeV) 


-325.1 


-179.3 


-128.1 


-123.0 


-357.7 


-250.3 


-296.7 


P (MeV) 


238.3 


71.9 


42.2 


51.6 


286.3 


149.6 


199.3 


7 


1.14 


1.35 


2.14 


2.14 


1.15 


1.19 


1.14 


-,„.(MeV fm^) 


21.8 


26.5 


18.3 


14.1 


19.5 


25.6 


22.9 


ir.(MeV fm^) 


-5.5 


-7.9 


-4.9 


-3.0 


-11.3 


0.0 


-2.7 


9r (MeV) 


5.9 


13.9 


6.4 


1.1 


0.0 


11.0 


9.9 


Csym (MeV) 


30.1 


33.0 


28.2 


27.0 


30.9 


29.3 


32.0 


asym (MeV) 


62.4 


29.8 


38.9 


42.9 


94.0 


79.3 


130.6 


hsy^ (MeV) 


-38.3 


-5.9 


-18.4 


-22.0 


-63.5 


-58.2 


-123.7 


Csym (MeV) 


-6.4 


-3.0 


-3.8 


-5.5 


-13.0 


-4.1 


12.8 


Poo (fm"^) 


0.16 


0.155 


0.145 


0.144 


0.162 


0.16 


0.16 


Jfoo (MeV) 


215 


262 


353 


366 


200 


239 


230 



mass of the pions. The Coulomb interaction is given by 

Nb 

VP^^^^Y.—^ (23) 

j = l J 

where the N-^ and iV^ are the total numbers of pions and 
baryons including charged resonances. The pion optical 
potential Rel/^^* originates from the medium effects in 
the hot and dense nuclear matter. In the calculation, we 
can also choose that of the vacuum, i.e., the Rey^J'P* is 
set equal zero. The influence of the pionic mean field 
in heavy-ion collisions on the transverse momentum dis- 
tribution was investigated by using a phenomenological 
ansatz and a microscopic approach based on the A-hole 
model by Fuchs et al. (20j . Here we use the phenomeno- 
logical ansatz suggested by Gale and Kapusta 21] . Then 
the dispersion relation reads 

u{p„ p,) = V^(|pJ-po)2+m2 - U, (24) 

U = ^pl + ml-m^, (25) 

mo = TOtt + 6.5(1 — a;"'")TO7r, (26) 
pI = {I- xfml + 2mam-^{l- x). (27) 

The phenomenological medium dependence on the 
baryon density is introduced via the coefficient x{pi) = 
exp(—a{pi/po)) with the parameter a = 0.154 and the 
saturation density po in nuclear matter. Influence of the 
in-medium effects on the charged pion ratio is also inves- 
tigated in Ref. 22]. 

III. RESULTS AND DISCUSSIONS 
A. Charge distributions 

Nuclear multifragmentation is a common phenomenon 
observed in heavy-ion collisions at intermediate energies, 
in which one can learn more about the properties of the 



fermionic nuclear matter, such as the compressibility, the 
liquid-gas phase transition and also the isospin asymmet- 
ric EoS. As a test of the ImlQMD model, we calculated 
the charge distributions in multifragmentation reactions 
of the system ^^^Au-|-^^^Au at different incident energies 
labeled in Fig. 2. The fragment multiplicity decreases 
with the atomic number and the trend is more rapidly 
at the higher incident energy The available experimental 
data [23, 24] can reproduced rather well besides the en- 
ergy at 35 A MeV in the region of the intermediate mass 
fragments (IMFs) (Z=3-15). The failure description of 
the production of the IMFs at several tens A MeV is due 
to the roughly treatment of fermionic nature in the dy- 
namical evolution of nucleons. These reactions are well 
described by the Statistical Multifragmentation Model 
(SMM) in which one assumes the formation of a thermal 
source and uses the statistical approach for its breakup 
[is! . Combination of transport models and SMM can re- 
produce the experimental data at fermionic energies. We 
constructed the fragments with the coalescence model, in 
which nucleons of the reaction system are considered to 
belong to a cluster in the phase space with the relative 
momentum smaller than Pq and with the relative dis- 
tance smaller than Rq (here Po=200 MeV/c and Ro=2A 
fm). 



B. Production of the total pion and the tt /tt^ yields 

Production of pion in heavy-ion collisions in the region 
of lA GeV is a primary product in nucleon-nucleon in- 
elastic scattering and mainly produced by the decay of 
the resonance A (1232). So it has enough yields to mea- 
sure pion in experiments to get the information of the 
EoS. Dynamics of the pion emission calculated by trans- 
port models is helpful for understanding the experimental 
observables. We calculated the time evolution of the mul- 
tiplicities of the total tt, A(1232), A^*(1440) and the ratio 
p/po in central ^^^Au-|-^^^Au collisions at lA GeV inci- 
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FIG. 2. (Color online) Comparison of the distribution of charged fragments in central ^^^Au+^^^Au collisions with the available 
experimental data at different incident energies. 



dent energy as shown in Fig. 3. One can see that pions 
are mainly in the domain at supra-saturation densities of 
compressed nuclear matter larger than the normal den- 
sity po • The production of pions is also influenced by the 
Fermi motion of baryons in the vicinity of the threshold 
energies. Figure 4 is a comparison of the measured total 
pion multiplicity by the FOPI collaboration in the reac- 
tion ''°Ca-|-*°Ca for head on collisions j4,] and the results 
calculated by the ImlQMD model for Skyrme parameters 
SkP, SLy6, Ska and SIII in the left panel, which corre- 
spond to different modulus of the incompressibility as 
listed in table 1, and the ratio of the calculated results to 
the experimental data at different incident energies. The 
total multiplicity of pion is mainly determined by the 
cross sections of the channels NN o A^A. The ImlQMD 
model with four Skyrme parameters predicts rather well 
the total yields at higher incident energies, but slightly 
overestimates the values near threshold energies, which 
may be influenced by the in-medium cross sections. In 
the calculation we use the in-vacuum cross sections of 
nucleon-nucleon elastic and inelastic collisions. Reason- 
able consideration of the in-medium inelastic collisions in 
producing A and N* is still an open problem in trans- 
port models, which have been performed in Giessen-BUU 
model 0. 

The 7T~ /7r+ ratio a sensitive probe to extract the high- 
density behavior of the symmetry energy per energy. Us- 
ing the isobar model 27|, one gets the ratio 7r^/7r+=2.09 
for pions from the A resonance, and 7r~/7r+=1.8 from 
the N* for the system ^'^^Sn-l-^^^Sn. These relations 
are determined by the Clebsch-Gordan coefficients and 
globally valid, i.e. independent of the pion energy. On 
the other hand, the statistical model predicts that the 
7r^/7r+ ratio is sensitive to the difference in the chem- 
ical potentials of neutrons and protons by the relation 
tt~/tt+ (X exp[2(^„-^p)/T] = ex.p[8SEsymip)/T], where 




10 20 30 40 50 60 

t (fm/c) 



FIG. 3. Time evolution of the total vr, A(1232), iV*(1440) 

ratio p/po in the reaction 
collisions at incident energy lA GeV. 



and the ratio p/po in the reaction ^^^Au+^^'^Au for central 



the T is the nuclear temperature [28]. The observed en- 
ergy dependence of the 7r^/7r+ ratio is due to the re- 
scattering and absorption process of pions and nucle- 
ons in the mean field of the compressed nuclear mat- 
ter. We use the free absorption cross sections in colli- 
sions of pions and nucleons by fitting the experimental 
data. The branch ratio of the charged tt and 7r° is deter- 
mined by the Clebsch-Gordan coefficients with the decay 
of the resonances A(1232) and N*(1440). The 77-/77+ 
ratio is sensitive to the stiffness of the symmetry energy 
at the lower incident energies. The compressed nuclear 
matter with central density about two times of the nor- 
mal density is formed in heavy-ion collisions in the 1 A 
GeV region. To extract more information of the sym- 
metry energy in heavy-ion collisions from the pion pro- 
duction, in Fig. 5 we calculated the 77-/77+ ratios as a 
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FIG. 4. (Color online) The excitation functions of the total pion multiplicities and the ratio of the calculated results with 
different Skyrme parameters to the experimental data in the reaction *''Ca+'"'Ca for head on collisions. 



function of the N/Z of the Sn isotopes in the reactions 
i00Sn+"2Sn, ii2Sn+ii2Sn, ^^^Sn+^^^Sn, i24Sn+i24Sn 
and ^^^Sn+^^^Sn in the left panel, and the excitation 
functions of the 7r^/7r+ ratios with the force SLy6 for 
the coefficients by Eqs. (12)-(15) besides the symme- 
try energy term. Here, the symmetry energy is treated 
as the cases without the term or with but different stiff- 
ness of the potential part which corresponds to superhard 
(7s=3), hard (7s=2), linear (73=!), soft (7s=0.5) and su- 
persoft (SHI). The symmetry energy reduces the 7r^/7r+ 
ratio for neutron-rich systems at the lower incident en- 
ergies because of its attractive force on protons and op- 
posite for neutrons in the dynamical evolution. So the 
probabilities of neutron-neutron collisions decrease after 
including the symmetry energy, which lead to the reduc- 
tion of TT^ production. The ratio increases with the stiff- 
ness of the symmetry energy, but decreases for the case 
of the superhard symmetry energy. The phenomena can 
be explained from the fact that although the symmetry 
energy reduces the N/Z ratio in the high-density region. 
A hard symmetry energy repulses neutrons of the collid- 
ing partners to flow the high-density region and enhances 
the domain of the supra-saturation density. Therefore, 
the neutron-neutron collisions in the case of the hard 
symmetry energy take place in a larger region at supra- 
saturation densities than a soft one. A too hard symme- 
try energy also reduces the value of tt^/tt"*" ratio because 
of the lower probabilities of the neutron-neutron colli- 
sions. The decrease of the Tr~ /tt'^ ratio with the incident 
energy is mainly owing to the production of pions from 
secondary nucleon-nucleon collisions, such as a neutron 
converts a proton by producing 7r~. Subsequent colli- 
sions of the energetic proton can convert again to neutron 
by producing tt"*". The 7r^/7r+ ratios of different situa- 
tion of the symmetry energy at high incident energies 
(Eiab >1.0A GeV) are close each other and approach the 



value calculated by using the isobar model. The larger 
■K~ /n^ ratio in the case of the hard symmetry energy 
was also reported in Refs 0, 01 ■ Shown in Fig. 6 is the 
profile of pion production without the symmetry energy 
and with increasing the stiffness of the potential part of 
the symmetry energy in Eq. (20). The symmetry energy 
reduces the tt" production owing to the decrease of the 
N/Z ratio in the high-density region. Influence of the 
stiffness of the symmetry energy on 7r+ is negligible. 

C. Influence of the in-medium pion dispersion 
relation 

Propagation of pion in the compressed hadronic matter 
is a complicated process that is involved in the interac- 
tion of pions and nucleons or pions and resonances. Fur- 
thermore, the pion in-medium effect has a non- negligible 
influence on the charged pion ratios. Calculations based 
on the QMD model show that the transverse momentum 
spectrum is not influenced with an inclusion of the A-hole 
potential [l^] ■ Opposite in the BUU calculations with the 
A-hole potential, the production of pion is enhanced at 
the low kinetic energies [i^. In Fig. 7 we show the influ- 
ence of the Coulomb potential and the symmetry energy 
in the mean field as well as the pion optical potential in 
the nuclear medium on the distributions of the transverse 
momenta of the charged and neutral pion production in 
the left and right panels, respectively, and also compared 
with the experimental data by FOPI collaboration [30| . 
We get the in-medium pion dispersion relation by using 
the phenomenological approach in Eqs (24)- (27). The 
propagation of pion in the nuclear matter is governed by 
the Coulomb interaction and the pion optical potential. 
Here we use the linear (7^=!) symmetry energy and the 
potential parameters in Eqs (12)-(15) are got from the 
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(N/Z) E, . (GeV/nucleon) 



FIG. 5. (Color online) The ■k~ j-n^ yields as functions of the neutron over proton N/Z of reaction systems in collisions of Sn 
isotopes at energy Eiai,=0.4A GeV and of the incident energy for the reaction ^'^^Sn+^^''Sn for the cases of different stiffness of 
the symmetry energy and without symmetry energy. 
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FIG. 6. Time evolution of the pion production for the cases of 
the linear, hard and superhard symmetry energy and without 
the symmetry energy in the reaction ^'^^Sn+^'^'*Sn at incident 
energy Eia6=0.4A GeV. 



Skyrme force SLy6 in table 1. The evolutions of produced 
pions in the baryonic matter are governed by the optical 
potential and also the Coulomb interaction of the charged 
pions and baryons. One can see that the Coulomb po- 
tential and the symmetry energy as shown in Eqs (9) and 
(20) have minor influence on the spectrum, especially 
for the TT"*" and 7r° production. The Coulomb interaction 
reduces the high-momentum tt^ production. The optical 
potential enhances the high-momentum pion production 
and also the tt"*" in the region of low comparing with 
the cases of the Coulomb potential and the symmetry 
energy. The same phenomena is observed in the QMD 
calculations for the tt'^ transverse momentum spectrum 



[ioj . The behavior can be understood by the strong at- 
tractive potential which enforces the pions to follow the 
trajectories of nucleons. Most pions get bound by the 
stopped participant matter resulting in an enhancement 
of the low -pt yield. On the other hand, pions which are 
bound by the spectator matter are driven out to high 
transverse momenta by the nucleonic flow and thus the 
high p( range is strongly enhanced. The different trends 
of the 7r~ and tt"*" are due to the Coulomb interaction of 
the charged pions and baryons. Calculations without the 
pion potential underpredict the high pion production 
is related to the missed resonances higher than N*(1440). 
The high-pt pion yields are described without introduc- 
ing pion potential by the calculations of Ref. 2^ which 
takes into account an extended set of baryonic resonances 
with masses below 2 GeV together with the in-mcdium 
reduced NN— J^NR cross sections. 

Figure 8 is a comparison of the longitudinal rapidity 
distribution with the same cases as shown in Fig. 7 for 
the reaction ^^''Au-|-^^''Au collisions at incident energy 
E/ab= lA GeV. A wider distribution is observed for the 
7r~ in the case of the optical potential. It is resulted from 
the fact that the optical potential enhances the high- 
energy pions especially for the 7r~ owing to the strongly 
attractive potentials. The optical potential underpredicts 
the mid-rapidity charged pion production. The symme- 
try energy reduces the mid-rapidity 7r~ production. The 
influence of the Coulomb potential and the symmetry en- 
ergy on the 7r+ spectra is slightly. We also calculated the 
total pion production and the 7r^/7r+ yields as a func- 
tion of the incident energies at the cases of the influence 
of the Coulomb, symmetry energy and pion optical po- 
tential, and also compared with the FOPI data Q as 
shown in Fig. 9. One can see that the total pion mul- 
tiplicity is not changed obviously by including any one 
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(MeV/c) p, (MeV/c) 



FIG. 7. Comparison of the transverse momentum distributions of charged pions and Tv° in the reaction ^^^Au+^^'^Au at incident 
energy Eiab= lA GeV for the cases of including aU, parts and none of the Coulomb, symmetry energy and the pion optical 
potential with the available experimental data [301 ]. 



or all parts of these terms. However, the 7r~/7r+ ratio 
is sensitive to the Coulomb potential and the symmetry 
energy owing to their sensitive dependence of the isospin 
degree of freedom of colliding systems, and weakly de- 
pends on the optical potential. The lower values of the 
7r^/7r+ ratios after including the in-medium pion disper- 
sion at near threshold energies are consistent with the 
analysis in Ref. 22|. Reasonable inclusion of the in- 



medium pion optical potential in the transport models is 
still an open problem, which is of importance not only in 
analyzing experimental data of the pion spectrum, also 
in constraining the high-density behavior of the nuclear 
symmetry energy. 



D. Flow distributions 

It is well known that the EoS at extreme conditions, 
the in-medium properties and the nuclear dynamics in 
heavy-ion collisions have been widely studied through 
the analysis of collective flow in experimentally associ- 
ated with transport models, such as nucleonic flow, light 
particles flow, and meson flow etc [1, [sil [s^] . It is pos- 
sible to reconstruct the reaction plane with flow analysis 
and hence to study azimuthal correlations. The analysis 
of the pion collective flow would be possible to get the in- 
formation of the dense isospin asymmetric nuclear matter 
formed in heavy-ion collisions. The directed flow vi is re- 
lated to the azimuth </> as wi = {Px/pt) = (cos^), and the 
elliptic flow f 2 = {{jPx/ptY-{py/ptY) = (cos(2(/))), where 
the angle brackets indicate averaging over all simulation 

pI- 



events and the transverse momentum pt = J pi: 



We compared the directed flows of the charged pions with 
protons and also with the FOPI data ^| as a function of 



the longitudinal rapidity in the ^^^Au-|-^^^Au reaction at 
incident energy 5A GeV as shown in Fig. 10. 

The rapidity is scaled by the incident projectile in the 
center of mass (cm.) frame as the relation = y/Up- 
The experimental data are selected in terms of the scaled 
impact parameter 6o = b/bma,x from the measured dif- 
ferential cross sections using a geometrical sharp-cut ap- 
proximation with &inax = 

1.15(Ap/^ + A^/^). We should 
note that the experimental data are sharply cut in the 
transverse momentum range 1.0< uto <4.2 scaled by pro- 
jectile four-velocity in the cm. frame Up as the relation 
Uto — Ut/up with the transverse four- velocity Ut- Calcu- 
lations are performed in the near central collisions with 
impact parameters b=2-4 fm and in the peripheral col- 
lisions with b=8-10 fm with the linear symmetry energy 
and without including the pion optical potential, but the 
Coulomb interaction of the charged pions and baryons 
is considered. The well-known 'S-shape' is clear in the 
distributions of the directed flows. The diagrams rotate 
clockwise with increasing the impact parameter. Calcu- 
lated distribution trends are consistent with the experi- 
mental results in both cases. Antiflow of tt"*" comparing 
the proton at larger impact parameter appears in both 
of the experiments and the calculations. The phenomena 
is caused by the shadowing effect of participant nucle- 
ons in heavy-ion collisions and was also investigated by 
BUU model [Hi and IQMD The difference of the 

flow of the charged pions is due to the Coulomb inter- 
actions. In Fig. 11 we present the transverse momen- 
tum distributions of the directed flow of charged pions 
and compared the recent FOPI data Q in the rapidity 
bin -1.8 < yo < 0. The ImlQMD model predicts the 
directed flow disappears with increasing the transverse 
momentum (uto >4). Calculated distribution trends are 
consistent with the experimental results in both cases, al- 
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FIG. 8. The same as in Fig. 7, but for the rapidity distributions. 




FIG. 9. (Color online) The excitation functions of the calculated pion multiplicity and the tt^/tt^ ratios in the reaction 
'^^'^Au+'^^'^Au for head on collisions with including all, parts and none of the Coulomb, symmetry energy and in-medium pion 
optical potential, and compared with the FOPI data i^. 



though the theoretical and experimental impact param- 
eter ranges do not precisely match. 

The elliptic flow in the ^^^Au-f ^^^Au reaction is also in- 
vestigated with the same framework of the directed flow 
as shown in Fig. 12. The flow spectra can be well re- 
produced in the near central collisions. However, the 
calculations overpredict the experimental data in the pe- 
ripheral collisions owing to the different definition of the 
impact parameters. The structures of the flow distribu- 
tions are similar. To explore more information of the 
isospin effect, we calculated the flow difference of the 
charged pious in Fig. 13, for that the isospin differ- 
ence is defined by Dvi = v'^ — v'^ for the directed 
fiow and Dv2 — wj^ — for the elliptic flow. One can 
see that the Coulomb potential plays a signiflcant role 
in the distribution of the flow difference. The influence 



of the symmetry energy on the distribution is slightly 
for the directed flow, and relatively apparent for the el- 
liptic fiow. But the symmetry energy at any stiffness 
does not change the distribution structure. The depen- 
dence of flow distributions of the charged pious on the 
Coulomb potential and on the stiffness of the symmetry 
energy was also investigated by the UrQMD model and 
the same conclusions were reported [35| . 



E. Polar anisotropy 

The analysis of the polar or azimuthal angle distribu- 
tions of pion emission in heavy-ion collisions is useful to 
extract the information of the hot and dense nuclear mat- 
ter. Anisotropic distribution of the polar angle of pion in 
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FIG. 10. (Color online) Comparison the directed flow of charged pions and protons with the FOPI data in the ^^^Au+^^^Au 
reaction at incident energy Eia6 = 1.5A GeV for the cases of near central and peripheral collisions. 
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FIG. 11. (Color online) Transverse momentum dependence of the directed flow of charged pions for the system ^^''Au-l-^^^Au 
at incident energy E(a(, = 1.5A GeV for the cases of near central and peripheral collisions, and compared with the FOPI data 
within the rapidity bin —1.8 < yo < 0. 



center of mass system was reported by a number of exper- 
iments [i, 36|. Isotropic pion emission is often assumed 
in the thermal model analysis of the compressed hadronic 
matter [33|. We calculated the polar angle distributions 
of the 7r~ in the reaction ^^^Au-f ^^^Au at different col- 
hsion centralities and also at the cases of the different 
stiffness of the symmetry energy as shown in Fig. 14. 
One can see that the anisotropy is pronounced in pe- 
ripheral heavy-ion collisions, which is consistent with the 
experimental observable [3|. The symmetry energies of 
the different stiffness have a minor influence on the po- 
lar angle distributions. Relatively, a wider distribution 
in the case of supersoft symmetry energy appears. The 
polar angle anisotropics result from the fact that pions 
are produced and emitted close to the surface of the re- 



action zone and they can escape without absorption by 
nuclear medium. An ellipsoid configuration in the re- 
action zone is formed in peripheral heavy-ion collisions 
and a significant fraction in the total pions are produced 
near the surface, so the large polar angle anisotropics are 
observed. 



IV. CONCLUSIONS 

The dynamics of the pion production in heavy-ion col- 
lisions in the region of lA GeV energies is investigated 
systematically by using the ImlQMD model. The calcu- 
lated total pions are consistent with the different Skyrme 
parameters SkP, SLy6, Ska, SIII which correspond to the 
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FIG. 12. (Color online) The same as in Fig. 10, but for the elliptic flow. 




FIG. 13. (Color online) Rapidity dependence of the flow difference of the tt^ and tt for the system ^^''Au+^^'^Au at incident 
energy Eiab=1.5A GeV and impact parameters b=8-10 fm. 
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FIG. 14. (Color online) Polar angle distributions of tt mesons in the reaction ^^^Au-|-^^^Au at incident energy Eiab~i-5A GeV 
for different collision centralities (left panel) and stiffness of the symmetry energy (right panel). 
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different modulus of incompressibility of symmetric nu- 
clear matter, and the available experimental data are re- 
produced rather well. The excitation functions of the 
■K~ /■n'^ ratio are sensitive to the stiffness of the symme- 
try energy. The difference is pronounced at near thresh- 
old energies and for neutron-rich systems. Calculations 
show that the ratio does not monotonically depend on 
the stiffness of the symmetry energy at supra-saturation 
densities. The influence of the in-medium pion dispersion 
relation on the distributions of transverse momentum and 
longitudinal rapidity is pronounced. But it has a minor 
contribution on the excitation functions of the total pion 
production and the Tr^/Tr"*" ratio. Calculations of the 
structure of the directed and elliptic flows are similar to 
the available experimental data. The flow distributions 
are sensitive to the Coulomb potential in the mean fleld, 
but weakly depend on the different cases of the symme- 
try energies. Polar angle anisotropics are pronounced in 
the peripheral collisions and the influence of the stiffness 
of the symmetry energy on the polar angle distribution 
is weak. 
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